Exact results for one-dimensional disordered bosons with strong repulsion 
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We study one-dimensional disordered bosons with strong repulsive interactions. A Bose-Fermi 
mapping expresses this problem in terms of non-interacting Anderson-localized fermions, whereby 
known results for the distribution function of the local density of states, the spectral statistics, and 
density-density correlations can be transferred to this new domain of applicability. We show that 
disorder destroys bosonic quasi-long-range order by calculating the momentum distribution, and 
comment on the experimental observability of these predictions in ultracold atomic gases. 

PACS numbers: 03.75.-b, 05.30.Jp 
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The properties of interacting bosons have recently at- 
tracted considerable attention due to the unprecedented 
control and tunability achieved in ultracold atomic gases. 
For instance, by using optical lattices, the predicted 
quantum phase transition from a superfluid state to a 
Mott insulator P, Q has been experimentally observed 
jsj. Current interest is also directed towards disordered 
systems, where disorder can be generated by using laser 
speckle patterns P|, additional incommensurate optical 
lattice potentials j^, iS] , or via atom-surface interactions 
in micro-chip confined atomic gases 0. It thus appears 
feasible to experimentally study dirty bosons under con- 
trolled conditions, in contrast to earlier realizations us- 
ing granular superconductors or Helium-4 in porous me- 
dia. Unfortunately, the theory of disordered interacting 
bosons is difficult, and no exact solutions are known apart 
from numer ical or approximate results ^.'^.'^IJinllTlIlT^. 
even in the one-dimensional (ID) limit .1?^ il4l |. 

In this paper, we show that for strong repulsive inter- 
actions, the dirty boson problem in ID is exactly solvable 
via a Bose-Fermi mapping discussed before in the clean 
limit of a Tonks-Girardeau gas [H [H [H [H [H 113 . 
The mapping establishes a connection to non-interacting 
disordered fermions, allowing to directly apply many re- 
sults previously obtained on Anderson localization in ID, 
and greatly simplifying the computation of other quan- 
tities like the momentum distribution. Our predictions 
can be checked using state-of-the-art experiments. De- 
tailed conditions for the ID regime have been specified 
in Refs. 0,122, and the ID Tonks-Girardeau regime has 
recently been achieved jlS,!^' ^1^° Refs. [25ll2^ . 

The Bose-Fermi mapping can be established most di- 
rectly by starting from a lattice description of hard-core 
bosons, the Bose- Hubbard model I'l, which applies im- 
mediately to optical-lattice experiments upon expanding 
the Bose field operator in the Wannier state basis 0. 
Considering spinless bosons on a ID lattice with spac- 
ing a, with Bose annihilation operator hi at site /, the 
Hamiltonian is 



blih 



em + —ni[ni 



where ni = bjbi- Here ei = hi + bP includes a random 
on-site energy hi and an axially confining harmonic po- 
tential, and Ji is a random hopping amplitude between 
neighboring sites. In optical lattices, hopping disorder 
is suppressed against on-site disorder j^, and we thus 
take Ji = J, but hi distributed according to a Gaussian 
ensemble |23 with 



0, hihi, = A6ii, 



(2) 



where the overbar denotes the disorder average and A the 
disorder strength. Using the spatial diffusion constant 



Dg, for a given disorder mechanism, A 



ti2„,3 



/(aDg) 



can be expressed in terms of microscopic parameters; vp 
is defined after Eq. 10 below. Detailed theoretical esti- 
mates for Dg (and hence A) are available for laser speckle 
fields and quasiperiodic optical lattices 0, , where 
also the Gaussian distribution of the disorder field hi is 
justified. We show below that the neglect of disorder in 
the Ji is no fundamental restriction. 

In the hard-core boson limit, U oo, only the oc- 
cupation numbers n; = or 1 are allowed, and then 
Eq. (0 can be mapped to a non-interacting lattice 
fermion model by means of a Jordan- Wigner transfor- 

- t 

mation, bi = e"^3<' '^j'^^ q, where the q denote lattice 
fermion operators. This transformation results in the 
fermionic Hamiltonian 



H 



I 



i+i 



^l^i Ci 



(3) 



It provides a one-to-one mapping, preserving the Hilbert 
space structure of the bosonic problem, with the A^- 
particle bosonic wavefunction expressed in terms of the 
fermionic one as |l6j| 



I In)\ 



(4) 
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The energy level Ei, for an A-boson eigenstate can 
thereby be computed in terms of the non-interacting 
fermionic Hamiltonian |(2Jl. In particular, with the single- 
particle energy for the jth fermion residing in a 
single-particle solution to Eq. (O, and taking into 
account the exclusion principle, = X^jLi 



The 



many-body fermionic wavefunction . . . ,In), and 

hence also the bosonic one (up to a sign), is then a 
Slater determinant, det['^i{lj)] /VnI. Since the mod- 
ulus square does not change under the mapping Q), all 
bosonic quantities given solely in terms of coincide 
with the fermionic ones. This includes all correlation 
functions of the particle density and the local density of 
states (LDoS), 

P(e'0 = E E S{e-E,)\<i>^{l,h,...jN)\^. 

V 1.2,. ...In 

The density of states (DoS) per site in a lattice with 
L sites is then p(t) — "^i^i pi^^l) / L, and also remains 
invariant. The same reasoning applies to the continuum 
limit studied later. In particular, the compressibility k, 
and thus also the sound velocity, is simply 

m \La / 

where m is the atomic mass. We set h — 1 and temper- 
ature to zero from now on. 

The equality of fermionic and bosonic results does not 
apply to the momentum distribution, 



II' 

Nevertheless, the Bose-Fermi mapping allows for a rather 
simple exact calculation of the disorder-averaged boson 
momentum distribution. Using the Jordan- Wigner trans- 
formation and Wick's theorem, Eq. (jSJ for a given dis- 
order realization can be written as a Toplitz determi- 
nant. For I > I', we find (6|5p) = 2'-''^Met[G('''')], 
where the (/ — I') x (I — /') matrix has the entries cfj ' = 
{cj,^^Ci,^j_i) — Sij-i/2, see also Ref. ^23|. For fixed disor- 
der {hi} and arbitrary trap potential, we compute Eq. ® 
numerically and subsequently average over different dis- 
order realizations. This is a much faster and more reliable 
procedure than directly studying interacting dirty bosons 
[l2.ll4| . since we have to deal with a single-particle prob- 
lem only. Let us first consider ®^Rb atoms in a harmonic 
axial trap with b — O.OIJ, see Eq. The overall energy 
scale is set by J, which can be tuned in optical lattices 
over a wide range (23| . We show results for N = 50 atoms 
in Figure n Clearly, disorder has a significant effect on 
the momentum distribution. In particular, some weight 
is transferred to large momenta, and the zero-momentum 
peak decreases, see inset in Fig. ^ The momentum dis- 
tribution has already been measured for bosonic atoms 
in the clean ID limit using Bragg spectroscopy [S^, and 
through imaging of the atom cloud after sudden removal 
of the trap potential |23| • Applying these techniques to 
the disordered case would allow to test our predictions. 
Note that these changes in the momentum distribution 
reflect Anderson localization physics, and should differ 
even qualitatively in the small-?/ limit. 
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FIG. 1: (Color online) Momentum distribution Q for several 
A and A'' = 50 rubidium atoms in a harmonic axial trap. For 
the disorder average, at most 300 disorder realizations were 
sufficient for convergence, and L was chosen large enough to 
ensure L-independence. Note the linear-logarithmic scale. In- 
set: Zero-momentum peak as a function of disorder strength. 



For a clean homogeneous (6 = 0) system, the bo- 
son momentum distribution is well-known to possess a 
h{p ^ 0) cx singularity 0,0], corresponding to 

the one-particle density matrix p{x,x') oc \x — x'\~^^^ for 
|a; — a;'| — s- oo. In the thermodynamic limit, Bose- Einstein 
condensation is absent, but there is quasi- long-range or- 
der characterized by the law. Remarkably, disorder 
has a fundamental effect on this singular behavior. The 
reasoning of Ref. 29] allows to prove that n(0) now must 
remain finite. The full momentum distribution n{p) is 
obtained numerically and shown for a ring with periodic 
boundary conditions in Fig. |21 The complete destruction 
of quasi-long-range order by disorder is clearly visible, 
and the momentum distribution becomes remarkably flat 
for sufficiently strong disorder. 

The Bose-Fermi mapping can also be established via 
the low-energy theory |13| , which is a perhaps more nat- 
ural description for magnetically trapped or micro-chip 
confined atoms, where no underlying lattice is present 
0. Focussing on circular or hard- wall axial trap poten- 
tials (6 = 0), corresponding to periodic or open bound- 
ary conditions, the resulting fermionic theory coincides 
with the continuum limit (a 0) of Eq. For a 

generic incommensurate filling N/L, we decompose the 
operator q into right- and left-moving {■ipR,'ipL) compo- 
nents with momenta k « ifep = ±7r iV/La according 
to Q ~ y/a[e^'^''^iljji{x) + e~'^''''^ipL{x)], where x = la. 
Correspondingly, the random on-site energies hi can be 
decomposed into a slow part and a term varying on 
a microscopic scale, hi « fj.{x) + {£,{x)e~^'^^''^ + h.c). 
For hi — 0, we have a ID massless Dirac Hamiltonian, 
f(x) produces a complex- valued random mass term, and 
li{x) a random chemical potential. Additional disor- 
der in the J; can be included in iji{x) and ^(cc), and 
with bispinor ip = (V'iJjV'i)j the continuum model reads 
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FIG. 2: (Color online) Momentum distribution on a ring for 
various A, N/L = 1/2, with L = 600. The A = result is 
analytical and shows the p scaling at p — » 0. Note 

the double-logarithmic scale. The inset shows that for finite 
A and L > 200, finite-size effects are negligible. 



H = J dx hip with 



(7) 



where vp = 2aJsin(fci?a) is the Fermi velocity. Here the 
2x2 matrices ct^ — [a^ ± icr^)/2 are defined in terms of 
Pauli matrices cr* acting in spinor space. Equation (UJ is 
the standard non-interacting Hamiltonian used to s tudy 
Anderson locahzation in ID conductors [23, ISfl IsJ ■ 
The forward scattering term can be eliminated by a 
gauge transformation and does not affect the quantities 
of interest below. (This is not possible on half-filling, 
where ^(a;) is real- valued and important differences arise 
(23.) Equation (0) then imphes 



f(x) = e(.^) = o, 



x'), (8) 



with T = i/vF for mean free path I, which also gives 
the localization length. We only discuss weak disorder, 
kpi :S> 1, where the bosonic system is in the Bose glass 
phase The average DoS is simply p{e) = \/{-kvf) 

|27| . and we turn to the LDoS probability distribution. 

The bosonic LDoS (from now on normalized to the 
average DoS) can be expressed in terms of eigenstates 
^i{x) with energy of the Hamiltonian Q, p{e,x) = 
TTVFj2i l*i(2;)P<5(£ — Cj)- ^ finite and closed sample, 
these levels are discrete and sharp, and it is necessary 
to regularize the J-functions. A natural way [3lj | is to 
smear out the S-peaks, pf{e,x) — J de' p{e' , x)f{e—e'), for 
instance by using a Lorentzian weight function, /,,(e) = 
r7/(7r[e^ -I- ry^]). Physically, the width -q is determined 
by inelastic processes, finite sample lifetimes, and escape 
rates of the trap. For an infinite sample, pf then follows 
the inverse Gaussian probability distribution |3]j| 



W{pf) 




g-4r,r(p/-l)Vp/ 



(9) 



which decreases exponentially both as a function of p/ 
for Pf 00, and as a function of 1/pf for pf — > 0. The 
anomalously small probability to find small pf implies a 
Poisson distribution of the energy levels, indicating the 
absence of correlations among close-by levels. This is 
obvious in the fermionic picture, where energy levels of 
localized non-overlapping states cannot repel each other, 
fn the strongly interacting bosonic picture, where well- 
defined single-particle states need not exist, this is a much 
less obvious result. Let us average pf{e,x) also over a 
spatial range S determined by the spatial resolution, e.g., 
the wavelength of a probe laser. For l/kp S <^ £ and 
ArjT ^ 1, the resulting LDoS p{e,x) is independent of 6 
and follows the distribution 31] 



Wip) 



rjT 



dt t smijTriTt) 



t + 4 
i-4 



rjrt 



(10) 

which is a somewhat narrower distribution than W{pf). 
Both Eqs. lO and (|10|l remain valid also in a finite closed 
sample, as long as the distance to any boundary is large 
compared to i. The LDoS can be measured usingirnaging 
methods or two-photon Bragg spectroscopy ^JEJl ; thus 
allowing for experimental checks. 

We then briefly discuss the bosonic LDoS correlations 
R{u;,x) at different energies and locations. 



R{lu, X — x') — /5(e, x)p{e + to, x') — 1, 



(11) 



which equal the fermionic ones computed in Ref. |3C| . 
The correlator pifl describes fluctuations in the spec- 
tral statistics related to energy level repulsion or at- 
traction. It is translationally invariant and independent 
of the energy e after the disorder average [s^- Since 
R{u!, X = x') we consider kplx — x'\ ^ 1 

and OJT ^ 1, where the correlator is finite, with the 
limiting values R = —1/3 for small x and i? = for 
X 00. In addition, there is a deep dip for £ <^ x <^ 
zo{u}) — 2£ln(8/a;T), where 



R{uj,£^ a; < zo) = -1 + 



16(x/£)3/2' 



(12) 



Here zo{uj) is the distance two nearly degenerate local- 
ized states must have to generate the energy splitting lu. 
The dip (|12|) implies that two states with nearly equal 
energies occupy with high probability locations far away 
from each other. Nevertheless, the wavefunctions of these 
states must have an appreciable overlap for short dis- 
tances, X < £. For X > zq, the LDoS correlations ap- 
proach the uncorrelated limit, 



R{ljJ, X > Zq) 



erf 



(13) 



where erf(x) denotes the error function. These features 
illustrate that the localized states are centered on many 
defects, leading to a complicated quantum interference 
phenomenon. As a consequence, close-lying levels do not 
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obey the usual Wigner-Dyson spectral correlations found 
in granular metals (s^l, but instead follow the Poisson 
statistics of uncorrelated energy levels. It would obvi- 
ously be quite exciting to probe experimentally. For 
ultracold bosonic atoms, this is possible usin g sti mulated 
two-photon Bragg scattering spectroscopy [23 • Simi- 
larly, one can compute the (Fourier-transformed) bosonic 
density-density correlations K(x — x' ,0;), 



Ke K{x,uj+iO) ~ IT J de nF{e-)[l—nF{e+Lo)] p{x; e+uj, e), 

(14) 

where npit) is the Fermi function, and the fermionic 
spectral function 

p{x ~ x';e + Ll>, e) — 

Y^S{e + uj- e,)6{e - e,)^,{x)^*{x')^*ix')^,ix) 



is a phase-sensitive quantity without a direct bosonic im- 
age. This nicely illustrates that the Bose-Fermi mapping 
opens otherwise unavailable routes for calculation. For 
X > zq [3^. 



while for x zo,to very good approximation [soj . p{x; e+ 
iu,e)^{TTVF)-^{R{uj,x) + l). 



To conclude, we have provided exact results for 
strongly repulsive dirty bosons in ID, which can be ob- 
tained from a Bose-Fermi mapping to non-interacting dis- 
ordered fermions. A Bose-glass phase is thereby mapped 
to an Anderson-localized fermionic phase. A similar map- 
ping is also available for arbitrary interaction strength, 
but involves interacting fermions with a non-standard 
contact interaction 0|. For strong (but finite) repulsive 
bosonic interactions, the weak fermionic interactions can 
safely be treated on a perturbative level jl^ and cause no 
substantial differences to our predictions. Finally, other 
quantities not discussed here can also be inferred, e.g., 
the crossover from short-time diffusive wave-packet ex- 
pansion to localized behavior at long times |35|. 



p{x; € + UJ,€) 



exp 
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